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ABSTRACT 

We investigate dynamical interactions between turbulent convection and g-mode pulsations in ZZ Ceti variables 
(DAVs). Since our understanding of turbulence is rudimentary, we are compelled to settle for order of magnitude 
results. A key feature of these interactions is that convective response times are much shorter than pulsation peri- 
ods. Thus the dynamical interactions enforce near uniform horizontal velocity inside the convection zone. They 
also give rise to a narrow shear layer in the region of convective overshoot at the top of the radiative interior. 
Turbulent damping inside the convection zone is negligible for all modes, but that in the region of convective 
overshoot may be significant for a few long period modes near the red edge of the instability strip. These con- 
clusions are in accord with those reached earlier by Brickhill. Our major new result concerns nonlinear damping 
arising from the Kelvin-Helmholtz instability of the aforementioned shear layer. Amplitudes of overstable modes 
saturate where dissipation due to this instability balances excitation by convective driving. This mechanism of 
amplitude saturation is most effective for long period modes, and it may play an important role in defining the red 
edge of the instability strip. 



1. INTRODUCTION 

The current paper focuses on dynamical interactions between 
turbulent convection and g-mode pulsations in ZZ Ceti vari- 
ables. It is the second of a series, and the last in which we 
rely on order of magnitude reasoning and the quasiadiabatic 



approximation. The reader is referred to Goldreich & Wu 1998 



(hereafter Paper I) for descriptions of the properties of g-modes 
and scaling relations appropriate to turbulent convection. Only 
results of immediate relevance to the current investigation are 
quoted here. Symbols taken from Paper I are defined in Table 
1. 

ZZ Ceti variables possess surface convection zones. Thus the 
interaction of their g-modes with turbulent convection requires 
attention. An important feature of convection in these stars is 
that its response time is much shorter than the periods of the 
observed modes. Thus to a good approximation, the convec- 
tion adjusts to the instantaneous pulsational state. This behavior 
plays an essentia l role in the overstability mechanism (Brickhill 
1990 , Gautschy 1996 , Paper I) which, following Brickhill, we 
refer to as convective driving. 

The nonlinear advection term in the fluid momentum equa- 
tion describes interactions between g-modes and turbulence. 
We model these interactions as giving rise to a turbulent viscos- 
ity which acts to reduce the magnitude of a mode's shear tensor. 
We devote §2 to evaluating the effects of turbulent viscosity on 
g-modes, both inside the convection zone and in the region of 
convective overshoot at the top of the radiative interior. We ap- 
ply these results in §3 to estimate the contributions of turbulent 
viscosity to mode damping. In §4, we analyze the saturation 
of the amplitude of an overstable mode due to instability of the 
narrow shear layer it forces at the top of the radiative interior. 
We conclude with a brief discussion in §5. 

2. EFFECTS OF TURBULENT VISCOSITY 



approximation. Perturbation theory is done with Lagrangian 
variables which are assumed to carry a time dependence e~ lM 
and a horizontal spatial dependence e' kh ' x . Gravity modes in- 
volve a mostly horizontal sloshing of fluid. Our principal con- 
cern is with the vertical derivative of the horizontal component 
of the displacement, d^h / dz, because it is the largest component 
of the grad ient of the displacement. 

In §2.1 we calculate the magnitude of d^hjdz inside the 
convection zone. To begin, we estimate the magnitude that 
d£h/dz would have in the absence of turbulent viscosity. In 
the dual limit of an isentropic convection zone and a vanish- 
ing entropy perturbation, \d£j l /dz\ would be equal to £\£ Z \/R 
which is smaller than £\^h\/R. However, a real convection zone 
is not an isentrope; its specific entropy increases with depth. 
Moreover, a g-mode produces a specific entropy perturbation. 
It turns out that the latter departure from isentropy is more im- 
portant than the former in enhancing the magnitude of d^h/dz. 
Without turbulent viscosity, the entropy perturbation could in- 
duce a \d^h/dz\ as large as £a/z w . Inclusion of turbulent viscos- 
ity re duce s \d£h/dz\ by a multiplicative factor cjf cv <C 1. 

In § 2.2 we examine the velocity shear in the transition region 
between the base of the convection zone and the top of the ra- 
diative interior. The turbulent viscosity declines sharply with 
depth in this region. We approximate this decline by a discon- 
tinuous drop to zero viscosity. The application of appropriate 
boundary conditions then shows that the horizontal velocity is 
discontinuous at Zb- We evaluate the magnitude of this velocity 
jump and the associated jump in the gradient of the pressure 
perturbation. 

2.1. Velocity Shear in the Convection Zone 

The linearized equations of mass and momentum conserva- 
tion, the latter in component form, read (cf. Paper I) 



$P ., c d£ z 

— = -IKhKh ~ —T- ) 

p dz 



Since motion in the outer layers of a DA white dwarf is all 
that concerns us, it proves convenient to adopt a plane-parallel 
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Table 1 
DEFINITIONS 



Symbol Meaning 

R stellar radius 

g surface gravity 

r radial distance from center of star 

z depth below photosphere 

Zb depth at bottom of convection zone 

Zuj depth at top of mode's cavity, z w <~ uj 2 /(gkl) 

uo radian mode frequency 

n radial order of mode 

I angular degree of mode 

kh horizontal wave vector, k\ = £(£ + 1 )/R 2 

p mass density 

p pressure 

s specific entropy in units of k B /m p 

F energy flux 

L luminosity, L = 4irR 2 F 

c s adiabatic sound speed, c 2 = (dp/dp) s 

p s (d\np/ds) p 

5 denotes Lagrangian perturbation 

£/, horizontal component of displacement vector 

£. vertical component of displacement vector 

v cv convective velocity, v cv <~ (F / ' p) 1 ^ 

f cv response time for convection, f cv ~ z/v cv 

v turbulent kinematic viscosity, v ~ zv cv 

A, B, C dimensionless constants approximately 2, 8 & 8 for ZZ Cetis 

Tth thermal constant at depth z, t cv /r t h ~ (v cv / c s ) 2 in the convection zone 

Tb unconventional thermal time constant at Zb, Tb ~ Tth/5 at Zb- 

t c time constant of low pass filter for convection zone, t c = (B + C)Tb 
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U & = m g£ z 



p p 

p d ( dp 
p dz 



Ja 
p 
dp 



+ ikh£,i 



A 
p : 



(2) 
(3) 



where /* is the force per unit volume due to turbulent convec- 
tion. 

For the moment, we neglect the force due to turbulence. Then 
differentiating equation (Eh and substituting for d£ z /dz using 
equation ([!]), we obtain 



dz w 



d (p Sp\ (5p 

-j- +8 — 

dz\P P J \ P 



With the aid of equation (|j) and the equation of state, 
Sp p Sp 



P 



-+p s 6s, 

CsP P 



we recast equation (|4j) in a more revealing form 

ighp. 



dz oj 



gpdz 



(4) 



(5) 



(6) 



We can think of the Sp/p and Ss terms as providing adiabatic 
and nonadiabatic forcing of d^/dz. Both terms vanish for adi- 
abatic perturbations in an isentropic convection zone. In their 
absence, the perturbations are irrotational so d^h/dz = ikh£, z - 3 
This value of shear would lead to negligible turbulent damping. 
To estimate d^h/dz when the flow is rotational, we must relate 
5p/ p and Ss to £/,. Here we appeal to Paper I which establishes 
that for z <Zb <z^, 4 



'111 
P 



and 



Ss ■ 



A(B + C) (bp 
1 — iu>Tc 



(7) 



The ratio of adiabatic to nonadiabatic forcing, as estimated at 
z = Zb, is given by 



p ds s 
8Pdz / 





dp 


1 


). 


p 


b \ 5s \b 



p_ds\ [l+(^r c ) 2 ]'/ 2 
gpdzj b A(B + C) 



+ Mcv) 



1/2 



«1, (9) 



where we make use of the relations gds/dz ~ (v cv /c s ) 2 and 
tcv/ T b ~ ( v cv/ c s) 2 - This establishes nonadiabatic forcing as the 
principal driver of velocity shear in the convection zone. In 
an inviscid convection zone, \dln^/dz\ would be of order z" 1 . 
Acting on a shear of this magnitude, turbulent damping would 
overwhelm convective driving. 5 

Next, we assess the effects of turbulence on the shear. Given 
our present understanding, about the best we can do is to model 
the turbulent Reynolds stress by analogy with the stress due to 

3 This is a linearized version of Kelvin's circulation theorem. 



molecular viscosity in a Newtonian fluid. The viscous force per 
unit volume is given by 



fi ■ 



daij_ 

dx-i ' 



(10) 



where the stress tensor, an, is related to the g-mode's displace- 
ment field by two scalar coefficients, the shear viscosity, v, and 
the bulk viscosity, according to 



-tujp < v 



dxi dx' 



3 at* u 



(11) 



When the effects of turbulent viscosity are neglected, d^/dz 
is by far the largest component of the gradient of the displace- 
ment. Since this term contributes to the shear but not to the 
divergence, we need not consider terms involving the bulk vis- 
cosity. Also, it is easy to show that f z /gp is much smaller than 
either Sp/p or Therefore, we ignore the effects of tur- 

bulent viscosity in the vertical component of the momentum 
equation. 

Turbulent viscosity has a profound effect on the horizontal 
component of the momentum equation. The dominant term in 
ft, takes the form 6 

f ~ • d ( d &\ ■ d & n ^ 

f h K.-ILU— pV— ~-lUipV CY — . (12) 

dz \ dz J dz 

Adding the approximate expression for —fh/p to the right hand 
side of equation (||) yields 



oj €h~ih g£ z + — — . 

p p J f C v dz 



(13) 



Differentiating this equation with respect to z, and repeating the 
steps that led to equation (0), we obtain 



(8) ^ 



LOtr- 



djh 

dz 1 — ioJt c 



gkhPs 



U3* 



gpdz 



(14) 



Since uit cv -C 1, turbulent viscosity significantly reduces 
d^i, I dz. 

For completeness, we provide an estimate for d£ z /dz- Com- 
bining equations ([[]), (§), (0), and (§) yields 



~ 77C/,^;, , 

dz 



(15) 



which is independent of the magnitude of the turbulent viscos- 
ity. The small value of d^Jdz ensures that its contribution to 
mode damping is negligible, as are the contributions from 
and ikh^z- 

2.2. Shear Layer at the Base of the Convection Zone 

Here we show that there is a shear layer at the top of the ra- 
diative interior across which the horizontal displacement jumps 
by an amount similar to the depth integrated change it would 
undergo in an inviscid convection zone. 



Overstable modes have Zlo > Zb- These are the only ones that concern us in this paper. 



A similar conclusion is arrived at by Brickhill ( 199C ). 



We neglect the depth variation of £/, with respect to that of pv. 



4 



We assume that the turbulent viscosity drops discontinuously 
to zero across the boundary at z = Zb- Then the continuity of the 
tangential stress demands that 7 



dz 



= 0. 



(16) 



b~ 



Next we multiply equation ( 13|) by p and integrate the resulting 
expression from the top of the stellar atmosphere to the bottom 
of the convection zone. This procedure yields 



dzp£,h 



' ikh I dzp 
Jo 



P\P 



(17) 



We transform equation (|l7]) in two steps. Integrating by parts, 
we obtain 



dzpih 



ik h pb 



dzp 



dz 



p 

p 



(18) 



Then using equations ([!]), (|3|), and (||), and taking into account 
the continuity of 5p and £, across Zb, we arrive at 8 



u> 



dzp£,h 



ui 2 p b 



61*+ + ih / dzp p 



p ds 
gpdz 



-Ss 
(19) 

Since for ojt cv <C 1, 6 is nearly constant within the convection 
zone, we establish the discontinuity of £/, across Zb to be 



-igki, 
uj 2 p b 



dz pp s 



p ds 
gpdz 



-Ss 



It follows directly from equation (pb that 



p d 
gpdz 



b* 



b- 



if 



(20) 



(21) 



Equation ( |20| ) merits a few comments. It may be simplified 
by taking advantage of the near constancy of Sp/p and Ss in- 
side the convection zone. Both are due to rapid mixing by tur- 
bulent convection which occurs on the time-scale f cv . Since the 
term involving the entropy perturbation dominates that due to 
the unperturbed entropy gradient, and \Ss\ ~ \Sp/p\ ~ 
the fractional change in £/, across Zb is of order Zb/zu- Thus the 
relative size of the jump is greatest for long period modes in 
cool stars. Finally, as advertised, comparison with equation (^) 
reveals that the discontinuity in £/, is closely related to the total 
variation that £/, would experience within the convection zone 
in the absence of turbulent viscosity. 

3. TURBULENT DAMPING 

Damping due to t urbu lent viscosity is evaluated inside the 
convectio n zo ne in §3.1, and in the region of convective over- 



shoot in §3.1 



3.1. Inside the Convection Zone 

Damping due to turbulent viscosity within the convection 
zone may be estimated as 



u?P? 



dzpv — — 
o V dz 



(22) 



Evaluating this integral with only the (dominant) entropy per- 
turbation term from equation (||) retained on the right hand side 
of equation dTi]), we obtain 



/vis- 



-0.1 



U> t rv T r 



1+(UT C ) 2 \Z 



Zb 



Sp 



(23) 



where (5p/p)b^ 1/ (rar w ) is the normalized eigenfunction eval- 
uated at Zb- In deriving equation (|23|), we integrate over an isen- 
tropic convection zone with adiabatic index 5/3, and substitute 
for A, B, C, v cv , and f cv according to Table 1. Radiative damping 
as evaluated in Paper I yields 



7rad 



-0.01L ( 6 -l 

V p 



(24) 



Since t cv /r c <C 0.1 for DA white dwarfs crossing the ZZ Ceti 
instability strip, turbulent damping inside the convection zone 
is much smaller than radiative damping for all overstable g- 
modes. 

It follows from equations (|l4|) and (^2|) that turbulent damp- 
ing of a given mode is maximized for ujt cv ~ 1 ; it is negligible 
inside the convection zone because uit cv <C 1 there. 10 However, 
ujtcv rises with depth in the region of convective overshoot. This 
suggests that damping in this region may exceed that inside the 
convection zone. 



3.2. In the Region Of Convective Overshoot 

The intensity of turbulence in this region diminishes with 
depth. We estimate its contribution to mode damping by solv- 
ing a model problem. We replace the bottom of the convection 
zone by a horizontal plate, and the region of convective over- 
shoot by an underlying viscous fluid of uniform density whose 
viscosity decays exponentially with depth, 



v(z) = v b e 



<z-z b )/X 



(25) 



The distance below the convection zone over which turbulent 
mixing maintains the entropy gradient small and negative is of 
order Aln(c s /v cv )|, or several times as large as A. 

The shear layer is driven by the plate which oscillates at fre- 
quency lu with displacement -A£/, given by (p0[). The Navier- 
Stokes equation yields 



d ( d^h 

TzVlTz ■ - l ^ = °> 



(26) 



The former is shown to be negligible for all 
modes. However, the latter may stabilize low frequency modes 
which would otherwise be overstable. 

7 We denote by feT quantities evaluated on the convective and radiative sides of the boundary. 
8 We drop the negligible lu 2 !;. and -/j/p terms in equation (^|). 

9 Here we treat the region of convective overshoot as having a finite extent unlike what we did in 

10 It might seem paradoxical, but provided tut cv <^ 1, turbulent damping is inversely proportional to the magnitude of the turbulent viscosity. This follows because 
the dissipation rate is proportional to (z 2 /f cv )(rf^;,/c/z) 2 , and \d^i,/dz\ oc LUt cv for uf cv <g 1. 



a 
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which is to be solved subject to the boundary conditions £/, = 
-A£/, at z = Zb, and £/, — > as z — > oo. Next we define 



«1, 



(27) 



(28) 



change independent variable to 

and dependent variable to 

w = e _fe-a)/2A f ft . (29) 
J) is recast in a more familiar 
1 



With these changes, equation 
form 

d 2 u 1 du 

— r+ + | 1 - 

dy z y dy 



r 



u = 0, 



(30) 



which is Bessel's equation for n = 1. The appropriate solu- 

(2) 

tion, which vanishes as \y\ — > oo, is ffj = 7i — 2Y1, where /i, 
Fi and i/{ 2) are the Bessel functions of the first kind, the second 
kind a nd the third kind, respectively (cf. Abramowitz & Stegun 
1970). The solution for £/, may be written as 



^ h =DyH^(y), 



(31) 



where the coefficient D is evaluated by applying the boundary 
condition £/, = -A£/, at z = Zb- Since \y b \ = 2e <C 1, we adopt the 



approximation ybH\ (y b ) ~ 2i/w, which implies 



D : 



inA£_ h 



(32) 



The viscous stress the plate exerts on the fluid is given by 



S = iujpbVb 



dz 



Making use of the identity 



i(2) 



(33) 



(34) 



we find 



~^X 



HhylH^(y b ) « 2cJ 2 \p b A£ h (lne + 7 + iV/4) : 



(35) 

where 7 is Euler's constant. 

The power per unit area supplied to the fluid, P/A, is equal 
to the product of S with the plate's velocity, iu>A^- To obtain 
the time averaged rate of energy dissipation, we multiply P/A 
by 4irR 2 and average over both time and solid angle to arrive at 



dE 
dt 



7TR 2 Pb L0 3 X\A^ 



(36) 



as it is in the interior of the convection zone. This rate first 
increases with depth, reaching a peak where the viscous diffu- 
sion time is of order and subsequently declines. Near the 
peak, v ~ ljX 2 , which implies a dissipation rate per unit vol- 
ume ~ pbOJ 3 \ A£/,| 2 . Taking the peak to have width ~ A, and 
including a factor R 2 for the integration over area, we recover 
equation (M) up to a factor of order unity. 

To relate^E' / dt) to the dampingrate due to convective over- 
shoot, we make use of equations (0), and (^|), and retain only 
the (dominant) entropy perturbation in equation (p0[). Adopting 
values and relations for A, B, C, f cv and v cv from Table 1, we 
arrive at 



-0.1- 



UJT C 



\ + {lot c ) 2 \z, 



A 



dp" 
P //• 



(37) 



Figure Q displays 7 v is-os/7rad for g-modes in stars of T e s = 
12,400K and T eS = 12,000K. n We choose \/z h = 1 in or- 
der to maximize turbulent damping in the region of convective 
overshoot. 12 Even with this extreme choice for X/zb, 7vis-os is 
smaller than 7 rac j for all the overstable g-modes. 

4. AMPLITUDE LIMITATION DUE TO SHEAR INSTABILITY 

Kelvin-Helmholtz instability of the shear layer at the top of 
the radiative interior is a source of nonlinear mode damping. 
As such, it acts to limit the amplitudes which overstable modes 
can achieve. We derive limiting amplitudes by balancing dissi- 
pation rates due to shear instability against excitation rates due 
to convective driving. In this section we deal with physical, as 
opposed to normalized, perturbations. 

Even in carefully controlled laboratory settings, unstable 
shear layers are complex structures. Those of interest here are 
further complicated by externally driven turbulence associated 
with penetrative convection. 13 In particular, turbulent mixing 
reduces the Brunt-Vaisala frequency while turbulent viscosity 
lowers the effective Reynolds number. The former reduces sta- 
bility while the latter enhances it. It is customary to express the 
nonlinear stress that maintains the shear in terms of the velocity 
jump across the layer, Av, and a di mensio nless dr ag coe fficient, 
Co, as (see, e.g., Landau & Lifshitz 1976, Tritton 1977) 



S= l -C D p{Av) 2 



(38) 



Co depends only weakly on Av in turbulent shear layers. Ter- 
restrial experiments indicate that Co varies logarithmically with 
the ratio of wall roughness to boundary layer width, with val- 
ues as small as 10~ 3 being characteristic of flows over smooth 
walls. We might speculate that penetrative convection makes 
the upper boundary of the shear layer behave like a rough wall. 
Since we have no physical basis for assigning a reliable value 
to Co, we treat it as a free parameter subject to the constraint 
10" 3 < Co < 10" 1 . 

The power per unit area input to the shear layer is P/A = 
SAv. To obtain the time average rate of energy dissipation in 
the shear layer, we multiply P/A by 4irR 2 , and carry out appro- 
priate averages over both time and solid angle. 14 Next we set 



Note that {dE /dt) is independent of v b and proportional to A. 
This is not surprising. Close to the plate the fluid moves almost 
rigidly, so the rate of dissipation per unit volume is small, just 

1 1 The DA white dwarf models used in this paper are provided by P. Bradley. For model details, see Bradley 
12 Helioseismology constrains the entire depth of the region of convective overshoot in the sun to be less than 0.05 pressure scale-heights (Basu 
13 We distinguish turbulence associated with convective overshoot from that due to instability of a g-mode's velocity shear. 
14 The velocity jump associated with a g-mode shear layer has an angular dependence and varies harmonically with time. 
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FIG. 1 . — Numerical solution to the toy model in § p.l| for e = 0.1. From top to bottom the panels display scaled dimensionless versions of the viscosity, v, the mag- 
nitude of the horizontal displacement, the magnitude of the principal component of shear, uj\d^i,/dz\, the local Reynolds number, Re = (ui\^i,\ 2 )/(u\d^i 1 /dz\), 
and the viscous rate of energy dissipation per unit mass, dE/dt = uui 2 \dS,i,/dz.\ 2 , as functions of distance below zt measured in units of A. Each of the curves plotted 
in the bottom four panels is a function of the single variable v/(io\ 2 ) and the parameter e. Moreover, the e dependence is weak (logarithmic) for e <C 1. Note that 
the shear, the Reynolds number, and the rate of energy dissipation all peak near the depth where*; /(u>\ 2 )= 1 . 
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FIG. 2. — Ratio of the rate of linear turbulent damping in the overshoot region to the rate of radiative damping for two white dwarf models with 7" c ff = 12, 400 K 
(solid triangles) and 12, 000 K (open squares), plotted against mode period for I = 1 modes. Here, A is taken to equal z.}, to maximize the effect of overshoot damping. 



Av = o;|A£/,|, where |A£/,| is evaluated using equation (0), 
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and (EOh as in §3.2. This series of steps leads to 



dE v 



dt 



OACd 



Zb L . (39) 



The saturation amplitude is obtained by balancing the damp- 
ing due to shear instability by the net rate of energy gain due to 
convective driving plus radiative damping (see Paper I), 



dE ne 
dt 



0.01 



[(^T C ) 2 -1] 
[K) 2 +l]' 



This yields 

Sp 
P 



o.i {^T C f+n l i\^T C ) 2 -nk h zi 



Cn 



LUT C 



Zb 



(40) 



(41) 



Amplitudes of overstable modes limited by shear instability 
exhibit the following patterns. In a given star, (5p/p)b declines 
sharply with increasing period. For a fixed mode, the value of 
{8p/p)b rises as the host star cools. The fractional flux per- 
turbation at the base of the convection zone is about 2(5 p/ p)b- 
The flux perturbation at the photosphere is smaller by the vis- 
ibility reduction factor 1/[1 + (wt c ) 2 ] 1 ' 2 (see Paper I). Realis- 
tic estimates of mode amplitudes must await consideration of 
additional amplitude limiting processes and the calculation of 
nonadiabatic growth rates. 

Is it safe to assume, as we have been doing, that the shear 
layer is turbulent? Necessary conditions are that the unper- 
turbed shear layer has both local shear large compared to u> 
and local Reynolds number large compared to unity. These 
conditions can be expressed by the dimensionless relations 
\d£ h /dz\ > 1 and Re = (w\£ h \ 2 )/(i/\d£ h /dz\) > 1. Based on the 



discussion in §3.2, we know that the left hand side of each in- 
equality attains a peak value ~ | A£& | / A near to where i/~wA 2 
(see Fig. []]). Thus the assumption of turbulence for low £ modes 
requires (5p/p)b ^> (_\/zb)(z^/R). The right hand side of this 
inequality is of order Zb/R ~ 10 -4 for lowest frequency modes 
detected in stars near the red edge of the instability strip. 

5. SUMMARY 



Brickhill (1990) discussed the manner in which turbulent 



convection affects g-modes in ZZ Ceti stars. He recognized that 
turbulent viscosity forces the horizontal velocity to be nearly 
independent of depth inside the convection zone. Moreover, 
he deduced that turbulent damping inside the convection zone 
is negligible whereas that in the region of convective overshoot 
might stabilize long period modes near the red edge of the insta- 
bility strip. Our investigation supports Brickhill's conclusions, 
although it suggests that turbulent damping by penetrative con- 
vection is, at best, of minor significance. 

We extend the investigation of turbulent dissipation to in- 
clude that due to instability in the shear layer at the top of the 
radiative interior. Turbulence generated by Kelvin-Helmholtz 
instability depends upon the velocity shear, so the rate of energy 
dissipation depends nonlinearly on the mode energy; dE /dt oc 
£ 3 / 2 . Thus shear instability provides an amplitude saturation 
mechanism for overstable modes. 

The applicability of results from this paper to ZZ Cetis must 
be viewed with caution. The quasiadiabatic approximation, 
while useful, is of limited validity. F ully n onadiabatic results 
will be presented in Wu & Goldreich (1998). 



We are indebted to Bradley for supplying us with models of 
DA white dwarfs. Financial support for this research was pro- 
vided by NSF grant 94-14232. 
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